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ABSTRACT 

We present a photometric study of the globular cluster (GC) system associated to the lenticular 
galaxy (SO) NGC 6861, which is located in a relatively low density environment. It is based 
on GEMINI/GMOS images in the filters g', r', i' of three fields, obtained under good see¬ 
ing conditions. Analyzing the colour-magnitude and colour-colour diagrams, we find a large 
number of GC candidates, which extends out to 100 kpc, and we estimate a total population of 
3000 ± 300 GCs. Besides the well known blue and red subpopulations, the colour distribution 
shows signs of the possible existence of a third subpopulation with intermediate colours. This 
could be interpreted as evidence of a past interaction or fusion event. Other signs of interac¬ 
tions presented by the galaxy, are the non-concentric isophotes and the asymmetric spatial dis¬ 
tribution of GC candidates with colours {g'- *')o>l .16. As observed in other galaxies, the red 
GCs show a steeper radial distribution than the blue GCs. In addition, the spatial distribution 
of these candidates exhibit strong signs of elongation. This feature is also detected in the inter¬ 
mediate subpopulation. On the other hand, the blue candidates show an excellent agreement 
with the X-ray surface brightness profile, outside 10 kpc. They also show a colour-luminosity 
relation (blue-tilt), similar to that observed in other galaxies. A new distance modulus has 
been estimated through the blue subpopulation, which is in good agreement with the previous 
value obtained through the surface brightness fluctuations method. The specific frequency of 
NGC 6861 {Sm = 10.6 ± 2.1) is probably one of the highest values obtained for an SO galaxy 
so far. 
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1 INTRODUCTION 


It is known that globular clusters (GCs) are usually ancient stel¬ 
lar systems, with age s that can be established within a reliable 
and tight scale (e.g.. ISalaris & Weisj|2002l : iMendel et al.|[2007ll 
that make them good tracers of the early evolutionary stages of 
their host galaxies. From the observational point of view, GC sys¬ 
tems reveal, to a greater or lesser extent, a bimodal c olour distri¬ 
bution reaching, in so me cases, a trimodal one (e.g. iBlom et d] 
l2012l:ICaso etalj|2013h . This multimodality indicates the presence 
of, at least, tw o subpopulations of GCs , usual l y refen'ed as ‘blue’ 
and ‘red’ (e.g. Gebhardt & Kissler-Patig|[l9^ : luarsen et~^l200ll : 


iKundu & Whitmore! 200 ih . They display different characteristics in 

terms of metallicity and projected spatial distribution, being th e red 
GCs more metallic jUsher et al.ll2012l : ISchuberth et alj|2010h and 
more concentrated towards their host galaxies than the blue ones. 


These differences suggest that both subpopulations arose 
through two different processes or in two phases. Among the most 
cited scenarios about the origin and evolution of GC systems and 


their h ost galaxies, we find: (i) the major merger formation sce¬ 
nario l lAshman & ZeDll[l992h which proposes that the merger of 
two spiral galaxies (S) leads to an elliptical galaxy (E) with ‘blue’ 
GC donated by the spirals and the ‘red’ clu sters formed in the 
merger; (ii) the multiphase collapse scenario l lForbes et al.lll99'J) 
which suggests that the ‘blue’ GCs are forme d first and t heir forma¬ 
tion is truncated by the cosmic reionization l lCerj|200lh . while the 
‘red’ GCs and the host galaxy field star s are formed late r, in a sec¬ 
ond phase; (iii) the accretion scenario jCote et al.lll998h . in which 
‘red’ GCs are formed in a massive seed galaxy and the ‘blue’ GCs 
are accreted from satellite galaxies of lower mass. T he inclusion of 
these scenarios in a proper cosmological framework ( Beaslevet^ 
I 2 OO 2 I; IPinino et alJlWl iMuratov & Gnedirj|201(]| ; lToninil|2013l) 
has shown different degrees of success and failure in explaining 
the observed properties of GC systems. Therefore, it is reasonable 
to assume that the formation process of GC systems may include 
ingredients from all of them. 

This close relationship between the galaxies and their associ¬ 
ated GC systems, allows us to compare the predictions of different 
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galaxy formation models with observational data. In this regard, 
we focus on the study of early-type galaxies and, in particular, 
of lenticular galaxies (SO) located in relatively low density envi¬ 
ronments. Most of the theories proposed for the origin of SOs are 
based on the re sults of different dynam ical processes on normal 
spiral galaxies dSnitzer & Baadel ll95lh. such as i nteraction with 
the intrac luster medium iCowi e & Songalij |l977h . ram-pressure 
stripping l Bekld^ al.ll2002l : Sun et alj 2006h . galactic harassment 
dMoore et a iJ I 19961) . etc. The presence of SO galaxies in low den¬ 
sity environments can hardly be explained just by these mecha¬ 
nisms. Therefore, they represent a challenge to our understanding 
of galaxy formation and evolution. 

In the last decade GCs were u sed to test the formation mecha¬ 
nism of SO galax ies in a few works dAragon-Salamanca et al.l2006l : 
iBarr et ani2007h . However, up to date, the number of lenticular 
galaxies with GC systems studied in detail remains relatively low. 
Therefore, in order to draw more firm conclusions, it is necessary 
to enlarge the sample of SO GC systems deeply analyzed. With 
that aims we started a study which is focused on several lenticu¬ 
lar galaxies belonging to the held or low density groups. 

In this context, we present here, as a hrst step, a study of 
the photometric properties of the GC system of the SO galaxy 
NGC6861 (Table Q}. This galaxy is located in the Teles copium 
Group (AS0851) at a distance of 28.1 Mpc dTonrv et alj|200ll) . 
and is one of the brightest galaxies in the group together with the 
E galaxy NGC6868. Chandra observations show that NGC6861 
might be undergoing an inter action with NGC6868 , suggesting 
the fusion of two sub-groups jMachacek et al.ll2010h . NGC6861 
presents an un usually high centra l stellar veloci ty dispersion of 
^414 km s“^ dWegner et alj|200^ . and from that iMachacek et ^ 
(l2010h inferred the pres ence of a supermassive black hole. Re¬ 
cently, iRusli et ^ ( l2013h included this galaxy in their sample of 
ten massive early type galaxies. They found that NGC 6861 (along 
with NGC 4751, another SO galaxy) shows the fastest rotation and 
the highest dispersion in their sample and confirmed the high dy¬ 
namical mass for this black hole (see Table[T]l. 

The GC system of N GC 6861 was previously studied by 
iKundu & Whitrnor^ ( 1200 ll) . They analyzed the GC systems of 29 
SO galaxies in the V and / bands (one image of 160 s in the F555W 
filter and two images of 160 s in the F814LP filter), observed be¬ 
tween 1995 September 6 and 1996 May 3 with the Wide Field and 
Planetary Camera 2 (WFPC2) of the Hubble Space Telescope. That 
work raises the possibility of the existence of bimodality in the 
colour distribution of the GC candidates in the range 0.5 <V-I< 
1.5 mag. Besides, the authors estimated a total GC population in 
their WFPC2 field of 1858 GCs and a corresponding local specific 
frequency of Sn = 3.6 ± 1.6. 



Figure 1. Digital Sky Survey image with the three superimposed GMOS 
fields. 


Table 2. Observations of NGC 6861, comparison field and blank sky. The 
values of airmass, exposure time and FWHM correspond to the final co¬ 
added images. Field 1 contains the galaxy. 


Galaxy 

Field 

Airmass 

g' r' i' 

9 

Texp. 

(sec) 

r 

i' 

FWHM 

(arcsec) 

9 r' i' 

NGC 6861 

1 

1.053 

1.056 

1.067 

4x700 

4x350 

4x400 

0.88 0.75 0.82 


2 

1.342 

1.251 

1.181 

4x700 

4x350 

4x400 

0.62 0.48 0.43 


3 

1.272 

1.095 

1.134 

4x700 

4x350 

4x400 

0.57 0.80 0.77 

Comp. Field 


1.026 

1.022 

1.039 

14x180 

14x120 

17x90 

0.56 0,58 0.53 

Blank sky 

- 

- 

- 

1.107 

- 

- 

7x300 

- - - 


those of the Sloan Digital Sky Survey (SDSS). The seeing condi¬ 
tions were excellent, ranging from 0.43 to 0.88 arcsec. As it can 
be seen in Figure[T] one of the fields was centred in NGC 6861. In 
order to fill in the gaps and facilitate removal of cosmic rays and 
bad pixels, four slightly dithered exposures were taken per filter. 

The raw data were processed by using GEMINI-GMOS rou¬ 
tines within IRAiQ (e.g. GPREPARE, GBIAS, GIELAT, GIREDUCE 
and GMOSAIC), and applying the appropriate bias and flat-field 
corrections. The bias and flat-field images were acquired from the 
Gemini Science Archive (GSA) as part of the standard GMOS 
baseline calibrations. 

It is known that i' and z' frames taken with GMOS show 
night-sky fringing. To subtract this pattern from our i' data, it 
was necessary to cieate fringe calibration images by using the task 
GIFRINGE. These frames were built from seven i' blank sky im¬ 
ages taken with exposure times of 300 seconds. Subsequently, the 
night-sky fringing was removed from our i' images using the task 
GIRMERINGE. 

As a final step, the reduced images corresponding to the same 
filter were co-added using the task IMCOADD, obtaining g', r' and 
i' co-added images for each field. These final images (see Table[3 
were then used for the subsequent photometric analysis. 


2 DATA 

2.1 Observations and Data Reduction 

We worked with images obtained in July-Nov ember 2010, using the 
Gemini Multi-Object Spectrograph (GMOS. [Hook et aill2004i of 
Gemini-South (Cerro Pachon, Chile) in imaging mode. The instru¬ 
ment consists of three 2048 x 4608 pixel CCDs, separated by two 
~2.8 arcsec (39 pixel) gaps, with a scale of 0.0727 arcsec pixel“^. 
The field of view (FOV) of the GMOS camera is ~ 5.5 x 5.5 
arcmin. Our images were taken under photometric conditions and 
using a 2 X 2 binning, which gives a scale of 0.146 arcsec pixel” ^ 
(Program GS-2010B-Q-2, PI: Bassino). Three de ep fiel ds were ob¬ 
served in the g', r' and i' filters dFukugita et al.|[l996l) . similar to 


2.2 Photometry and classification of unresolved sources 

The source detection was performed on the g' images because they 
present a better signal-to-noise ratio (S/N) than the r' and i' frames. 
In order to do that, we used a scr ipt that combines features of SEx- 
tractor dBertin & Arnoutslll996l) . along with IRAE filtering tasks. 
It allows us to model and subtract both the background and the 

^ IRAF is distributed by the National Optical Astronomical Observato¬ 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation 
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Table 1. Galaxy propierties. (2-5) Equatorial and galactic coordinates (NED) ; (6) morphological classification (NED); (7) B magnitude from the RC3 catalogue 
ide Vaucouleurs et alj|l99lh ; (8) distance m odulus from|Tonrv et alj i2QQll) ; (9) effective radius from RC3; (10) central velocity dis persion calculated into a 
circular aperture of radius r=0.595 h“^kpc jWegner et alj2003l) ; (11) velocity dispersion within Re and (12) black hole mass from jRusli et ^ .120131) . 


Galaxy 

“J2000 

5j2000 

1 

b 

Type 

Bt 

(m-M)o 

Re 

(^0 


Mbh 





(“ 


(mag) 

(mag) 

(") 

(km 

(km 

(Mq) 

NGC 6861 

20 : 07 : 19.5 

-48 : 22 : 13 

350 : 52 : 38 

-32 : 12 : 39 

SA(s)0“ 

11.92 

32.24 

17.7 

414±17 

388.8±2.6 

2.2(2.1,2.4)xl0® 


brightness due to the galaxy halo. To do this, initially, stellar ob¬ 
jects were removed by using SExtractor, and then the residual im¬ 
age was smoothed with a median filter, employing the FMEDIAN 
task from IRAF. This median filtered image was then subtracted 
from the original image. The procedure was applied a second time 
on the last image using a smaller median filter to detec t and to dis¬ 
card w eak halos around objects near the galactic centre dpaifer et al.l 
l201lh . For each field, the script provides a final catalogue with all 
detected objects and a galaxy light subtracted image. 

From these lists of objects, ~50 bright unresolved sources 
were selected per field, uniformly distributed over the FOV. Then, 
second order models of the point spr ead function (P SF) were ob¬ 
tained using the DAOPHOT package dStetsonll 19871) within IRAF. 
Several tests were performed with different psf models, selecting 
the moffat25 model, which gives lower fit errors than the Gaussian 
and moffatlS options. The same procedure was applied to the r' 
and i' images of the same field. In order to obtain the instrumental 
psf magnitudes of all the detected sources, we applied the ALLSTAR 
task on each image, using their respective psf models and the object 
list of the reference filter {g'). In addition, we conducted an aper¬ 
ture correction to the psf magnitudes through the MKAPFILE task. 

As a final step, we built a master photometric catalogue in which 
the coordinates of objects belonging to different fields were in the 
same reference system. To do this, we selected common sources 
located in the overlapping regions of the fields and worked with the 
tasks GEOMAP and GEOXYTRAN within IRAE. 

At the distance of NGC6861, we expect that GC candidates 
are unresolved sources. Therefore, the Stellarity Index of SExtrac¬ 
tor (0 for resolved objects and 1 for unresolved ones) was used to 
make the object classification. We set resolved/unresolved bound¬ 
ary in 0.5. 

2.3 Photometric Calibration 

Several standard star fields observed during the same nights as our 
frames were downloaded from the GSA data base (see Table O. 
All of them are from the list of ‘Southern Standard Stars for the 
u'g'r'i'z' System’ of Smith et al. (in preparation^ Exposure times 
for these images were 5.5 sec for the g', r' and 4.5 sec for the i' 
filter. They were reduced using the same bias and flat-field images 
previously used for the reduction (see Section IZTt . 

The PHOT task within IRAF was used to obtain instrumental 
magnitudes of the standard stars through large apertures of 3 
arcsec of radius. Thus, no aperture corrections were applied to this 
photometry. Subsequently, and with the aim of bringing our instru¬ 
mental magnitudes to the AB standard system, we used the trans¬ 
formation: 

tttstd — rttzero “i” C’T'jmi P^2')std (1) 

where rUstd are the standard magnitudes, mzero is the photometric 

^ http://www-star.fnal.gov 


zero point, rriinst are instrumental magnitudes, Kcp is the mean 
atmo^heric extinction at Cerro Pachon given by the Gemini web 
pagqj X is the airmass, CT is the colour term, and (mi — m2)atd 
is one of the co lours listed in the fourth column of Table|4l 

Following ljdrgens^ ) |2009[) . we applied the extinction correc¬ 
tion and got a mean magnitude zero point of zero order for each 
standard frame, filter and night (values listed in Tabled- After ob¬ 
taining these zero points, we applied them to our science photome¬ 
try and to the instrumental standard star magnitudes. Then, we used 
the full sample of these first order calibrated standard star magni¬ 
tudes to get the residual zero point for every star: 

XlTlzero — PP’std [mtnst Xqp^X 1) -f TTLzero] (2) 

We fitted these ‘residual’ zero points by using the following lin¬ 
ear relation: Arrizsro = a -|- CT{mi — m2)std- Figureshows 
Amzero versus the standard colours and the fitting functions. As 
can be seen in Table [T] in the case of g' and r' we got significant 
values of a and CT. However, the i' parameters are consistent with 
zero. Then, we applied the g' and r' colour terms to our previously 
calibrated science photometry. 

Finally, objects in common in each science field were used to 
obtain any small zero point differences between our field 2 and the 
other pointings, resulting values lower than 0.03 mag in all cases. 
These small offsets were applied in order to refer our photometry to 
field 2. Subsequently, we visually rejected spurious detections and 
artifacts from our master catalogue (less than 1 per cent of total 
sources) and applied the Gala ctic extinction coefficients, given by 
ISchlaflv & Finkbeined ( 1201 ih (A„ = 0.179, = 0.124, Ai = 

0.092). 

2.4 Completeness test 

To quantify a reliable detection limit of unresolved objects in our 
photometry, we performed a series of completeness tests in each 
filter and each field. We split the range 20.1 <go< 29.5 mag into 
intervals of 0.2 mag, adding 200 point sources using the ADDSTAR 
task. A total number of 9600 artificial objects were added in each 
field. The spatial distribution of the artificial objects was carried out 
using a power law, trying to follow the radial distribution of the GC 
candidates. Subsequently, the same script mentioned in Section l2(^ 
was applied to these images, in order to recover the added objects. 

Figure shows the completeness curves for the three fields, 
plotting the fraction of artificial objects recovered versus the input 
magnitude. It can be seen that, at go=25.5 mag, the data have a 
completeness greater than 80 per cent. Thus, in the following anal¬ 
ysis, we will take a faint cut magnitude of g'o = 25.5 mag. Figure|4] 
shows the completeness fraction for different ranges of magnitudes 
as a function of galactocentric radius (Rgai). It can be seen that 
most of the objects are lost at small Rgai, due to the high bright¬ 
ness of the galaxy. Taking this completeness variation into account, 
we can conclude that for magnitudes brighter than (;q= 25.5, we 

® http ://www. gemini. edu/?q=node/10445 
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Figure 2. Residual zero points as a function of the standard colour index. Red, grey and black circles represent the photometric measurements performed on 
standard stars in the fields mentioned in Table[3] The blue line shows the obtained fit. 


Table 3. Standard fields used for the photometric calibrations. *S: number 
of standard stars in the field. Kcp: atmospheric extinction coefficient at 


Cerro Pachon. 

zero point obtained from the fit. 

Standard field 

*s 

Science field Filter 

KcP 

’^zero 

020000.300600 

2 

1 9 ' 

0.18 

28.42±0.010 

220100300000 FI 

4 




180000.600000 

15 

r 

0.10 

28.33±0.003 

180000.600000 F5 

16 




E8.aFl 

12 




180000.600000 

15 

i' 

0.08 

27.92±0.004 

180000.600000 F5 

16 




E8.aFl 

12 




180000.600000 

15 

2 9 ' 

0.18 

28.41±0.006 

180000.600000 F5 

16 




220000395900 FI 

4 




180000.600000 

15 

r' 

0.10 

28.44±0.004 

180000.600000 F5 

16 




220000395900 FI 

4 




180000.600000 

15 

i' 

0.08 

28.01±0.006 

180000.600000 F5 

16 




220000395900 FI 

4 




020020.600000 

4 

3 9 ' 

0.18 

28.24±0.016 

CDFS 

2 




020020.600000 

4 

r 

0.10 

28.30±0.011 

CDFS 

2 




020020.600000 

4 

i' 

0.08 

27.92±0.005 

CDFS 

2 





Table 4. Residual zero point, colour terms and the final photometric zero 
points for field 2, obtained with the linear fits. 


Filter 

a 

CT 

(mi-m2) 

Mg ero 

g’ 

-0.054±0.019 

0.08ib0.03 

ia'-r') 

28.36±0.02 

r' 

-0.018±0.011 

0.03±0.01 

(g'-r') 

28.42±0.01 


0.006±0.013 

-0.02±0.05 

(r'-i') 

28.01±0.01 


detect the vast majority of GCs over the whole NGC6861 field, 
except in the inner 30 arcsec of Rgai (1-7 i?e) where most of the 
objects are lost. 



Figure 3. Completeness fraction as a function of Jq magnitude. The solid, 
dashed and dotted lines con'espond to fields 1, 2 and 3, respectively. The 
horizontal dotted line shows the 80 per cent completeness level. 


2.5 Comparison Field 

We were not able to get a comparison field as part of our Gem¬ 
ini program and, as we will mention in Section lTsI our GC candi¬ 
dates seem to fill the whole field. Therefore, in order to estimate the 
background contamination, we had to search for deep g', r' and i' 
fields in the GSA. We found only one nearby GMOS frame which 
filled our requirements. This field, centred in the compact galaxy 
group NGC6845, (RA=20'‘00™-58.28“, DEC=-47‘*04”‘11.9'’) is 
part of the Gemini program GS-201 lA-Q-81 (PI: Gimeno). We 
downloaded the raw images and reduced and calibrated them ac¬ 
cording to the procedure explained in Subsections l2.1l to l2.4l 

Unresolved objects were selected using the same ranges in 
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Figure 4. Completeness fraction as a function of galactocentric radius 
for field 1. The curves (from top to bottom) show the completeness for 
5 q= 22, 24, 25, 25.5, 26.5, 27 mag. The horizontal dashed line shows the 
80 per cent completeness level. 
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Figure 5. Colour-magnitude diagram of the points sources detected in the 
NGC 6845 field. 


colour and magnitude of our GC candidates in NGC 6861 (see Sec¬ 
tion [3]2j. In Figure]^ we show the colour-magnitude diagram for 
the unresolved sources obtained for this field. The completeness 
test on the g' frame shows that we have a completeness >80 per 
cent for g' < 26 mag. After that value, the completeness falls down 
strongly and colour errors grow rapidly. 

It is important to mention that the NGC 6845 group is com¬ 
posed of four galaxies. All of them are well contained inside the 
GMOS FOV. Two of these galaxies are spirals and two are of SO 
type. The only available estimation of the distance to the group is 
base d on the galaxy reds hifts and they give values from 90 to 100 
Mpc ([Gordon et al.l200^ . Adopting thi s range of distances and the 
V to g' band relation from iFaifer et al.l ( 1201 ih . we can expect that 
massive GCs like lo Cent will show g^ ~ 24.8 — 25.1 mag. Then, 
as a first step we should check if the GCs belonging to these galax¬ 
ies are detected in our photometry. The colour-magnitude diagram 



Figure 6. Colour-magnitude diagram of the points sources in the compari¬ 
son field. Black filled circles show colours and magnitudes inside the ranges 
adopted for the GC candidates in NGC 6861. 

shown in Figure|^ indicates that no particular grouping is detected 
for objects in the colour ranges typical of GCs. We performed the 
same test using the colour-colour diagrams (CDs) and the results 
were similar. As an additional test we made plots of the spatial dis¬ 
tribution of the point sources and only those objects with go >24.5 
mag seem to show some very marginal concentration towards the 
galaxies. Therefore, as expected from the assumed distance, some 
of the fainter objects in NGC 6845 may well be GCs belonging to 
the galaxies in the group. 

We excluded the zones near to the four galaxies and we ap¬ 
plied the same magnitude cuts as in the NGC 6861 sample (see Sec¬ 
tion and thus only 36 point sources survive in the NGC 6845 
field (Figurel^. This is a very low number and if we accept that 
some of them could be GCs belonging to the group, it means that 
our NGC 6861 photometry has a contamination lower than 8 per 
cent. In what follows we will take these 36 objects as our contami¬ 
nation estimation. 


3 RESULTS 

3.1 Galaxy surface profile and tidal structures 

We obtained the surface brightness profile of NGC 6861 using the 
IRAF task ELLIPSE which fits elliptical isophotes to the galaxy. In 
order to avoid the light contribution from nearby bright or extended 
objects, before executing ELLIPSE we masked them in the image. 
We modeled the galaxy light allowing the centre, ellipticity and po¬ 
sition angle {PA) of the isophotes to vary freely until the fit became 
unstable in the outer regions. It is important to remember that the 
GMOS FOV is small and NGC 6861 fills it completely. Therefore 
all the parameters mentioned here are limited to Rgai < 2 arcmin. 

The results are shown in Figure]?] It can be seen that the el¬ 
lipticity undergoes strong variations from 0.46 to 0.2 within 100 
arcsec, whereas the PA shows a difference of ~44° between 
the innermost zone and the radius mentioned above. This ellip¬ 
ticity decrement is also seen in the X-ray hot ga s emission from 
e=0.4±0.04 to e=0.19±0.05 for 7.3^r<42 kpc dMachacek et al.l 
[ 201 ^ . 

These changes in the isophotal parameters are also detected 





















6 Escudero et al. 



equivalent radius (arcmin) 


equivalent radius (arcmin) 


Figure 7. Isophotal parameters for the g' filter versus equivalent galactocentric radius. The latter was calculated as the square root of the product of the 
semi-major and semi-minor axis lengths of the fitted ellipses. From left to right and from top to bottom, the graphics show the ellipticity and position angle of 
the isophotes, and the Fourier coefficients A3, S3, A4, S4, respectively. 



Figure 8. Left: r' Isophotes of NGC6861 (field 1). The faintest and more external isophote correspond to /iy,/=24.5 mag/arcsec^, and brightest and inner to 
= .5 mag/arcsec^. The steps between them, are 0.368 mag/arcsec^. Right: Ratio between the the original r' image and the model. It is possible to see 

the different structures present inside of NGC6861. The most obvious and easily visible structures were labelled by letters ‘A’ to ‘E’. North is up and east is 
to the left. 


in the Fourier terms (A 3 ,S 3 ,A 4 ,S 4 ), which measure the isophote’s 
deviations from perfect ellipticity. The terms A 3 and give “egg- 
shaped” or “heart-sha ped” isophotes, but t he most significant term 
is the coefficient ,64 jjedrzeie^ 3 ll987l) . According to its sign, 
E/SO galaxies can be classified in disky (B^yO) or boxy (_B 4 < 0 ). 
In our case, within 1 arcmin. Bn indicates disky isophotes. At larger 
radii, the parameter B4 changes its sign showing that NGC6861 
displays boxy isophotes. These results are in agreement with those 


obtained bv iLi et ^ ( |20 1 ih in the The Camegie-Irvine Galaxy Sur¬ 
vey (CGS). It is interesting to men tion that hoxy isophotes are usu¬ 
ally suggestive of recent merging dKormendv & BendeJI 19961) . 

The va riation shown by the PA in NGC 6861 was previously 
reported by iTal et alj ( l2009h . These authors mentioned that this 
galaxy has tidal characteristics, showing intra-group emission and 
non spherical isophotes. Therefore we searched this kind of evi¬ 
dence in our deep GMOS images. In the left hand panel of Figure[ 8 ] 
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Figure 9. Radial surface brightness profile of NGC6861 in filter g'. 

we show the r' image of NGC 6861. It can be seen that the galaxy 
isophotes are not concentric, showing a deformation towards the 
N-W. The right hand panel in that figure shows the ratio between 
the original r' image and the smooth ellipse model. Several distinct 
structures are evident in NGC 6861, starting by a dust and stellar 
disk labelled with ‘A’ (Rgai<W arcsec), a series of more external 
arcs (where the brightest are labelled with ‘B’ and ‘C’), and two 
diffuse stellar debris or low surface brightness structures. The fea¬ 
ture located in the NW was labelled with ‘D’ and that in the SE 
with ‘E’. Unfortunately, ‘D’ presents a bright star on it. Between 
‘C’ and ‘D’ there is clearly another structure. The structure ‘E’, 
which is not symmetrically located with respect to ‘D’, is easily 
visible and shows a sharp external edge. All of them lie along the 
inner major axis of the galaxy. 

Einally, and in order to compare the spatial light distribution 
of NGC 6861 with that of the GCs, we obtained a more extended 
surface brightness profile for the galaxy. Although the images of 
fields 1 and 2 were taken on different days, Table|4]shows that the 
photometric conditions were similar. Therefore, we followed the 
same procedures as those included in the task IMCOADD, used the 
overlapping area of fields 1 and 2 to equalize the signal and the sky 
level of the g' images, and built a mosaic including these two fields. 
After checking the photometric scales by comparing common ob¬ 
jects in the original images with those in the final mosaic, we run 
ELLIPSE and we get a profile reaching 6 arcmin of sma (Eigure[^. 
The inner zone of the profile (sma<70 arcsec), which is saturated 
in our long exposure images, was obtained by running ELLIPSE on 
a short exposure g' image. 

3.2 Globular cluster colours 

FigurefTOl shows the g'o versus {g'-i')o colour-magnitude diagram 
(CMD), as well as the ig'-r')o versus {g'-i')o and the {g'-r')o 
versus (r'-r')o colour-colour diagrams, which include all the un¬ 
resolved sources detected in the three GMOS fields. As it can be 
seen, GCs are easily identified as point sources clustered around 
specific colours and most of these unresolved objects define a clear 
sequence, which is highlighted above the sequence of the Milky 
Way stars. In order to obtain our sample of GC candidates, we se¬ 
lected the objects within the colour ranges 0.5 <{g'-r')o< 0.95 



Figure 11. Colour-magnitude diagram of the final GC candidates (black 
filled circles), the dashed line shows the 80 per cent completeness level. 

mag; 0.05 <{r'-i')o< 0.5 mag; 0.55 <{g'-i')o< 1.45 mag. We 
expect that a small percentage of the GC candidat es are foreground 
stars and possible Ultra Compact Dwarfs (UCDsi lPrinkwater et al.l 
l2003h . To separate bright GCs from UCDs and/or field stars, we 
took Pq= 21.8 mag as the bright end in the CMD, which is roughly 
equivalen t to My —11 m ag, (in agreement with the value sug¬ 
gested bv iMieske et alJl^OOfih . The adopted faint end was gQ=25.5 
mag, where we ensure a completeness level greater than 80 per 
cent, and colour photometric errors lower than 0.1 mag. 

The final sample of GC candidates for the subsequent analy¬ 
sis is shown in Figure[TT] A total of 1245 objects meet the above 
criteria. The figure shows all the candidates identified before per¬ 
forming any statistical subtraction of field contamination. 

At first glance, the CMD seems to show bimodality, display¬ 
ing a marked blue subpopulation and a more disperse, but large, 
red subpopulation. Interestingly, FigurefTT] shows that the blue se¬ 
quence of GCs seems to become redder with increasing luminos¬ 
ity. This particular characteristic shown by NGC 6861, is known as 
‘blue-tilt’. We will discuss it in Subsection l3.4l 

3.3 Colour histograms 

Due to the large number of GC candidates, we were able to build 
colour histograms counting objects in small bins of 0.05 mag. In the 
upper panel of FigurefT^ we show the histogram obtained for all 
candidates in the sample and the expected contamination obtained, 
as it was indicated in Subsection l2.5l A factor of 3.0 was applied 
to the contaminant object counts in order to take into account the 
different areas of the NGC 6861 and NGC 6845 fields. As it was 
mentioned in that section, the number of MW stars and background 
galaxies in the GCs sample seems to be very low. Therefore, in the 
lower panel of FigurefT^we show the histogram corrected by the 
background contamination and, as expected, the correction had a 
negligible effect on the integrated colour distribution. In the same 
figures we show the respective smoothed colour distributions ob¬ 
tained using a Gaussian kernel of a—0.04 mag. This a value was 
considered to be representative of the mean error in {g'-i')o colours 
for fainter candidates. 

To quantify the modes and dispersions of the different sub¬ 
populations of GC candidates, we used the statistical algorithm 
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Figure 10. Colour-magnitude (a) and colour-colour diagrams (b,c) of all of the detected point sources in the three fields (black dots), and the initial sample 
of GC candidates (black filled circles). Mean photometric errors {g'-i')o are shown in panel (a) by red bar's. The red box shows the limits used in colour and 
magnitude for the selection of GC candidates. 


Gaussian Mixture Model (GMM: [Muratov & Gnedinir201(lh . This 
method is a parametric probability density function represented 
as a weighted sum of Gaussian component densities. Using three 
statistics, which are 1) parametric bootstrap method (low values 
indicate a multi-modal distribution), 2) separation of the peaks be¬ 
tween Gaussians (D>2 implies multi-modal distribution) and 3) 
kurtosis of the input distribution (k<0 condition necessary but not 
sufficient for bimodality), GMM quantifies whether a multi-modal 
distribution provides a better fit than a unimodal one. It also allows 
to define initial values for the modes and fits equal dispersions (ho- 
moscedastic) or different dispersions (heteroscedastic) to the com¬ 
ponents. 

We started by fitting two Gaussian distributions and finding 
mean colours for the expected ‘classic’ blue and red GCs. The 
brightest objects ((jio<22.8 mag) were rejected in order to avoid 
the effect of the ‘blue tilt’ and the usual unimodal colour distribu¬ 
tion of these massive objects. The values obtained, listed in Tabled 
are in good agreement w i th the location of the peaks found in other 
galaxies by iFaifer et^ ( 1201 ll) . In that work, the authors present 
results for another bulge dominated SO galaxy, NGC3115, which 
is considered now as one of the most clear examples of bimodal¬ 
ity in colour ((5'-i')o~0.765±0.007 and (g-t')o~ 1.044±0.011 
for bl ue and red GCs, respectively) and in metallicity dBrodie et al.l 
l2012h . Interestingly, NGC 3115 shows more or less the same lumi¬ 
nosity as NGC 6861, but the latter presents a larger GC population 
and its bimodality is not as clear. 

Figure[T2 shows that the red subpopula tion in NGC 6861 
presents a broad and flat colour distribution. iBlom et ^ l l2012h 
found similar characteristics in the massive elliptical galaxy 
NGC 4365, where the existence of a third subpopulation of GCs 
with intermediate colours is proposed. In order to study whether 
the same situation occurs in NGC 6861, we split the sample of GC 
candidates in four radial bins with approximately 310 objects in 
each of them (Rgai ^ 46, 46-93, 93-160 and >160 arcsec). The 
background-corrected colour distributions are shown in FigurefT^ 
A Gaussian kernel was used in order to obtain smoothed colour dis¬ 
tributions, as described earlier. We can see that in the inner region 
there is a clear blue subpopulation at modal value {g'-i')o ~ 0.85. 
Besides this, a clearer red peak is easily detected at modal colour 
of {g'-i')o ~ 1.10. As we move towards outer regions, and as 
expected, the blue GCs begin to dominate. However, the red GC 
candidates become more heterogeneous and at least two other pos¬ 


sible peaks seem to appear. One of them in {g'-i')o ~ 0.94 and the 
other one around {g'-i')o ~ 1-10 — 1.15 mag. 

With the purpose of quantifying and characterizing the possi¬ 
ble presence of three subpopulations, we ran GMM to the colour 
distribution and the RMDijj software to the histogram (whole sam¬ 
ple background-corrected histogram), fitting in both cases a tri- 
modal distribution. The obtained values for the GMM fit, listed in 
Table[3 indicate that a trimodal distribution is as good as a bimodal 
one. On the other hand, the reduced chi-square values obtained with 
RMIX indicate that a trimodal case gives a better fit to the data than 
the bimodal one. 

In different early-type galaxies, the mean colour of the red 
GCs were fou nd to be similar to t he inner halo colour of their 
host galaxies dForbes & FortelEoOlh . This is the consequence of 
the bul ge/spheroid field stars and red GCs having a strong genetic 
nexus ( Forte et ^l2009h . Through the task ELLIPSE we obtained 
an approximate colour profile for NGC 6861 (see Section lUTt . This 
profile was calculated within a radius less than 2 arcmin, in order 
to avoid the strong effect of the error in the sky level. In the up¬ 
per panel of Figure[T4] we compare the colour of the galaxy halo 
with the mean colour of the GCs for the bimodal case. To separate 
blue and red subpopulations, we consider the colour cut in {g'- 
i')o = 0.90. It can be seen that the red GCs look similar to the inner 
halo, although somewhat bluer. However, considering a trimodal 
case {{g'-i')o ^0.88, 0.88<(()'-T)o<1.01, and {g'-i')o ^1.01, 
for blue, green and red subpopulations, respectively), we can see 
that the red GCs look almost identical in colour to the halo (lower 
panel in Figure M- In both cases, the GC subpopulations were 
split using the values provided by GMM, and subsequently mean 
colours were obtained within several galactocentric intervals con¬ 
taining the same number of objects, 85 GCs in the bimodal case 
and 51 in the trim odal one. It is interesting to see that figure 16 in 
iBlom etal.l ( l2Q12h presents similar results. 

3.4 Blue tilt 

As mentioned in Section[T2] the blue GC candidates of NGC 6861 
seem to follow a colour-magnitude trend, known as ‘blue-tilt’. A 
similar phenomenon was observed in several elliptical galaxies 

RMIX is publicly available at http://www.math.mcmaster.ca/peter/mix/mix.html 
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Table 5. Values obtained with GMM and RMIX for the colour distribution. Statistics values (p, D, k, x^) listed along with the mean colour p, 

their dispersion cr and the fraction of objects assigned by RMIX for each subpopulation for the bimodal and trimodal case. 


Population 
Statistics values 


GMM 

Bimodal case Trimodal case 

(p=0.01; D=2.63±0.14; fc=-0.934) (p=0.01; D=2.46±0.34; A:=-0.934) 

li <7 fl (T 


RMIX 

Bimodal case IVimodal case 

(X^=I5.639) (x^=5.869) 

IJ, a f II a- f 


Blue 0.79±0.02 0.08±0.007 0.783±0.014 0.075±0.009 0.783±0.01 0.07±0.006 0.376±0.023 0.778±0.02 0.07±0.003 0.379±0.025 

Green - - 0.954±0.062 0.066±0.032 - - - 0.955±0.03 0.07±0.Q20 0.227±0.140 

Red 1.07±0.02 0.12±0.007 1.128±0.058 0.10±0.022 1.066±0.014 0.13±0.008 0.624±0.023 1.136±0.04 0.10±0.016 0.394±0.127 



Figure 12. Upper panel: colour histograms for the GC candidates 

and comparison field (dashed lines) with magnitudes between 21.8 <Sq< 
25.5. The solid and dotted lines represent the smoothed colour distributions. 
Bottom panel: Background corrected colour histogram and smoothed dis¬ 
tribution. 


(see Forte et al j2007l : IWehner et all2008l : lHarrisl200^ : |Peng et al.l 


I 2 OO 8 I) . This trend has been interpreted as an increase of the 
metallicity with the GC mass, i.e., a mass-metallicity relation 
(MMR). Physics explan atio ns for the ‘blue- t ilt’ we re proposed by 
I Strader & SmithI 12(3081) and iBailin & Harris l l2009l) . based on self¬ 
enrichment during a cluster’s formation stage, where star formation 
is influenced by supernova feedback within the protocluster. 

With the aim to characterize this MMR in NGC 6861, the sam¬ 
ple of GC candidates was split in several magnitude bins of 100 



(9'-i’)o 



(g’-Oo 


Figure 13. Colour distribution of the GC candidates, separated in four radial 
bins, according to their galactocentric distance {0-46, 46-93, 93-160 and 
>160 arcsec). The solid line represents the applied Gaussian kernel. 
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Figure 14. Colour index {g'-i')o vs. projected galactocentric radius. The 
bimodal case is plotted in the upper panel and the trimodal case in the bot¬ 
tom panel. Red circles, blue squares and green triangles represents mean 
values in colour for the red, blue and intermediate subpopulations, respec¬ 
tively. The black line show the profile colour of the galaxy within a 2.5 
aremin radius. 
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objects. Subsequently, RMIX was used to obtain the locus of peaks 
in each bin for the three subpopulations mentioned in section 1331 
(GMM gave us similar results, but RMIX was preferred due to its 
graphical direct interface). In this case, we made homosedastic 
runs in order to facilitate the convergence (the three modes have 
the same variance) and tested with bins of different wides. In Fig- 
ure[T5]we show, as an example, the mean values for blue, green and 
red candidates obtained considering bins of 0.075 mag in an i' ver¬ 
sus {g'-i')o CMD. The figure shows that, as was suspected earlier 
from the colour-magnitude diagram, the blue peak indeed becomes 
redd er at higher lum inosities. 

iHarrij ( l2009ah found evidence that the blue-tilt in M87 begins 
to be detectable from magnitudes brighter than Mi=— 9.5, exhibit¬ 
ing a n onlinear behaviour. Using the equation (1) from iFaifer et al.l 
(I 2 OI ih and the adopted distance to NGC6861, we obtained the 
equivalent valu e i'n ^ 23.24 mag. Regarding the slope of the 
tilt obtained by iHarrisl it was d{g — i)/di = —0.021 ± 0.004. 
However, values sligh tly larger than that were obtained for other 
galaxies. For example. iFaifer et"^ ( 1201 ih found d{g' — i'^jdi' = 
—0.032 ± 0.003 (fo r luminosity M i^—9.5) for the giant ellipti¬ 
cal galaxy M60, and IWehner et al.l (l2008h found d{g' — i')/di' = 
-0.044 ± 0.011 (for luminosity M/sJ-O) for NGC 3311. 

In order to compare with the different values from the litera¬ 
ture, we made a linear least-squares fit to the colour as a function 
of g'o and i'o magnitudes. We use this functional form as a first or¬ 
der approach and because it is worthless to use a more complex 
functional form for this particular sample of points. 

As an example, we list here the fit obtained considering all 
the measured points, that obtained rejecting the points fainter than 
i'o = 23.3 mag (Mi ~ —9.5 mag), and that obtained rejecting only 
the brightest point, respectively: 

{g' - i')o = -0.026(±0.007)io + 1.385(±0.163) (3) 

{g - i')o = -0.025(±0.011)io -b 1.358(±0.268) (4) 

{g - 1)0 = -0.025(±0.009)io -b 1.344(±0.232) (5) 

Although the behaviour of the data at the faint end is not so clear, 

we can see that the slopes for the different fits are identical. There¬ 
fore, we will adopt as the slope of the ‘blue tilt’ in NGC 6861, the 
first value, i. e., d{g' — i')/di'=—0.026(±0.007). 

In order to facilitate the comparison of our results with 
other values present in the literature, we additionally run 
RMIX forcing it to fit a himodai distribution to the same sam¬ 
ple used in the trimodai case. In this opportunity, we made 
heteroscedastic runs (that means that different mode variances 
are allowed) and the resulting peaks of the two Gaussians for 
different magnitude bins are also shown in Figure [HI We can 
see that the blue peak is clearly recovered, and that most of 
them are very close to the points obtained in the trimodai case. 
On the other hand, the new red peaks are located between 
the “green” and “red” candidates of the tri-modai case. A lin¬ 
ear fit to the blue peak points using the whole sample gives 
{g' - 1)0 = -0.021(±0.011)io -b 1.258(±0.257). When only 
the points brighter than i'o = 23.3 are included, the fit results 
in (g' - i')o = -0.033(±0.010)io -b 1.539(±0.239). Within the 
errors, these values are the same as those in equations 
and very similar to those obtained for several elliptical galax¬ 
ies, as mentioned before. 

The adopted values for the blue-tilt in the himodai or tri¬ 
modai cases, do not necessary imply that the MMR is linear for 



Figure 15. Blue-tilt present in NGC 6861 (black line) for trimodai case. 
Blue, green and red filled circles represent peak valnes in (g'-i')o 
colour obtained by RMIX for blue, green and red subpopniations, re¬ 
spectively. The brown open squares correspond to the peak valnes con¬ 
sidering the bimodal case for blue and red subpopulations. 


NGC 6861. Our tests indicate that the slope, considering different 
sub-samples, seems to decrease as we consider the weakest mag¬ 
nitude points in the CMD (as f ound in other gala xies like those 
presented by iMieske et aLlEoT^ and lHatTisll2009al) . However, the 
errors associated with those estimations are too big, indicating that 
we probably need a more deeper sample to obtain conclusive re¬ 
sults. For that reason, in this paper we decided to compare our re¬ 
sults only with estimates of linear fits in the literature. 

As mentioned before, the different estimates found in the lit¬ 
erature on the blue-tilt, mostly for elliptical galaxies, are similar to 
that obtained here. This is true even when we compare slopes from 
ot her combinat i ons o f photometric filters. For examp le, the results 
of Usher et ^ 1 20]J^for NGC 4278, or those from iMieske et al.l 
( l2006h and Istrader et alj ( l2006h for M49, M60 and M87, are all 
around d{g' — z')/dz^—O.G'iG. This value is equivalent to that ob¬ 
tained by us in NGC 6861 when w e use the re l ation between the 
ACS (g'-z') and {g'-i') filters from iForte et alj(l2013h . 

In a similar way as it was done for blue GCs, we per¬ 
formed linear fits for green and red candidates in the trimodai 
case, obtaining d{g' — i')/di'=—0.009{±0.022) and d{g' — 
i')/di' =—0.019(±0.018), respectively. As it can be seen, the 
slope for green GCs is not significant whereas in the case of the 
red GCs, a weak red-tilt is obtained, though the value of this slope 
has a marginal significance. When we force the bimodal case and 
we lit the read peak, we get d{g' — i')ldi' = 0.012(±0.009). 
This slope has the opposite sign to that of the trimodai case. 
However, this marginally significant slope disappear when the 
brightest point is rejected. Therefore, we think that no slope is 
present in the red peak sample when bimodality is assnmed. 


3.5 Spatial distribution. 

In Figure[T6lwe show the projected spatial distribution of the three 
subpopulations according to the separation of GMM mentioned in 
Section|33] It is clear in this figure the different spatial distributions 
of blue and red candidates. The latter looks more concentrated and 
flattened, resembling the starlight distribution of NGC 6861 (see 
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Table 6. Slope values for power law and de Vaucouleurs law, in the density 
profiles for the whole sample and for the blue, green and red subpopulations. 


Population 

Power law 

de Vaucouleurs 

All 

-!.60±0.04 

-0.87i0.03 

Blue 

-1.47±0.07 

-0.73i0.04 

Green 

-1.35±0.08 

-0.71 ±0.05 

Red 

-2.04±0.11 

-1.07±0.08 


Subsection 13. 6T On the other hand, we have obtained the radial 
unidimensional distributions of GC candidates. The surface density 
for each radial bin was obtained counting the GC candidates in an¬ 
nuli, and applying a factor to take into account the fraction of area 
effectively observed of each annulus. In every case, the associated 
density error is given by Poisson statistics. 

For the entire sample of candidates, we used concentric cir¬ 
cular annuli with A log r = 0.1, corrected by background contam¬ 
ination, reaching ~ 10 arcmin (~100 kpc) from the galactic cen¬ 
tre. Subsequently, we have obtained separated profiles for the three 
subpopulations: red, green and blue (Figure IT/I ). Due to the elon¬ 
gation showed by the red and green subpopulations, the profiles 
associated with them were obtained using concentric elliptical an¬ 
nuli, according to the position angle and ellipticity values obtained 
in Section [T6] In the case of blue GC candidates, we didn’t de¬ 
tect any significant elongation, therefore we decided to obtain their 
density profile using circular annuli. It can be seen in Figure [TT] 
that the red subpopulation presents a stronger spatial concentration 
towards the galaxy than the blue GCs. On the other hand, blue and 
green subpopulations show very similar radial profiles. Regarding 
the extension of the GC system of NGC6861, Figure [T^ shows 
candidates at galactocentric distances even larger than 10 arcmin 
(~100 kpc). Although at such distance our areal completeness is 
small, we can not rule out that some of them are bona fide GCs. 

In order to assess the slope of the profile, we fit a power law 
and a de Vaucouleurs law to the whole sample, and to the 

red, green, and blue candidates, respectively. We present the ob- 
tained values in Ta ble[^ These values are similar to those found in 
iFaifer et alj l l201lh for NGC 524 and NGC 3115, other SO galaxies 
(see their Table 5). It is interesting to note that the green candi¬ 
dates show identical slopes than the blue ones. However, Figure [771 
shows that green candidates display a possible excess of objects at 
the inner region of the system, not present in the blue or in the red 
profiles. 

In summary, from FigureslT^andl 17lwe can see that the spatial 
distribution for the blue clusters is radially shallow and circularly 
symmetric; for the green candidates is shallow and elliptical; and 
the red GCs distribution is stepper and elliptical. 

In order to see if there is any indication of sub-structure in the 
spatial distribution of our GC candidates, we analyzed the appear¬ 
ance of the xy plots by taking the candidates in different {g'-i') 
colour windows. That exercise has shown that GC candidates with 
{g'-i')>l.l6 present an asymmetric distribution around de centre 
ofNGC 6861. In the lower panel of Figurell6l we plotted those can¬ 
didates as black points. In order to quantify this asymmetry we have 
done an easy test; we have taken objects with 30 < Rgai < 160 
arcsec (complete areal coverage), and then counted them splitting 
the sample between the left (SE) and the right hand side of the 
galaxy (NW). We have got 69 and 35 candidates, respectively. 
Those numbers show that there are two times more such candi¬ 
dates to the SE of the galaxy centre than to the NW. This could 
be not a contamination effect because, according to the analysis 
of our comparison field, we can expect few contaminant object in 



Figure 16. Projected spatial distribution of GC candidates. From top to bot¬ 
tom, the blue, intermediate and red subpopulations, according to the GMM 
separation. Black filled circles are the GC candidates with (g'-i')>1.16 
which show an asymmetric distribution around the centre of NGC 6861. 
The orientation is the same as Figure[T] 



(R)^™ (arcsec) 

Figure 17. Projected density profile for red (triangles), green (diamond) 
and blue (squares) subpopulations, and all GC candidates (circles). Solid 
grey line shows the galaxy starlight profile, which has a similar slope to 
red clusters. All profiles were fitted with a power-law and a de Vaucouleurs 
law. The vertical dotted lines indicate the range used to make the fits (2.4< 
r-l/4 < 4.6 arcsec). The profiles were shifted to avoid overlapping. 


that area (<10). Inhomogeneities in GCs spatial distributions were 
reported in othe r early type galaxies lik e NGC 4261, NGC 4649 
and NGC 4278 (iD’Abrusco et al.ll2013alH l2014h . These features, 
like fossil remnants, are strongly suggestive of the accretion and/or 
merger of some lower mass neighbors. 
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Table 7. Values of position angle and ellipticity for the three subpopulations 
and for the whole sample of GC ca ndidates. In add i tion, w e include the 
results for X-ray measurements from iMachacek et 


Param. 

Blue 

Green 

Red 

All 

X 

r ^7.3 kpc r ^42 kpc 

PA 

e 


145°±13° 

0.54±0.1 

146.8°±5° 

0.64±0.04 

137°±6° 

0.52±0.05 

128°±7° 

0.40±0.04 0.19±0.05 


3.6 Azimuthal properties 

In Subsection 13.51 Figure [7^ we saw that red candidates show a 
flattened projected s patial distribution. In ord er to check that, we 
use the expression of lMcLaughlin et ^ (Il994h : 

a{R,9) = kR~°‘[cos^{9-PA) + {l-e^)~‘^sin{e-PA)]~°‘^'^ 

( 6 ) 

to study the GC azimuthal distribution, where cr is the number of 
GC candidates. We consider the whole sample, and the three sub¬ 
populations of GC candidates. We fit the position angle {PA), mea¬ 
sured counterclockwise from the north, and the ellipticity (e) in 
Equation]^ There, a is the value of the power-law exponent previ¬ 
ously obtained in the surface density fit (see Section [J31 >. In order 
to bring together the largest number of GC candidates, and avoiding 
correction for areal incompleteness, we took a circular ring around 
the galaxy with 31 <Rgai< 101 arcsec. Subsequently, we divided 
the ring in wedges of 22.5°, and we fit the Equation|6l allowing the 
normalization constant (k), the PA and the ellipticity to vary. In 
Figure rr^ we show the histograms of the azimuthal distributions 
for the whole sample, for the blue, green and red GC candidates 
and the obtained fits. In Tabled we present the values obtained 
for the two fitte d parameters. There w e also list the results for the 
X-ray emission l lMachacek et al.|[2010l) . 

Figure [T^shows that the spatial distribution of the whole sam¬ 
ple of GC candidates, exhibits strong signs of elongation. The PA 
results identical within the error to that of the galaxy inside a semi¬ 
major axis of 120 arcsec (see Subsec tion l3. lb . As far as ellipticity 
is concerned, the values obtained for green, red and all GC candi¬ 
dates are slightly bigger than that of the galaxy. In particular, red 
candidates show a significantly higher ellipticity than that of the 
starlight. It is interesting to notice that, according to Figure [T^ the 
behaviour of the GC candidates is led by the red and green ones. 
For blue candidates we didn’t find any significant elongation and, 
then, the PA results unconstrained. 

3.7 Comparison with the hot gas emission 

In their st udy of the GC sy stems of the giant elliptical galaxy 
NGC 1399. lForte et alj l l2005h have shown that there is a similarity 
between the projected density profile of the blue GC subpopula¬ 
tion and that o f the X-ray surface brightness profile of that galaxy. 
More recently. [Forbes et aP l l2012ll have extended the study to other 
eight giant elliptical galaxies and they have found that the same be¬ 
haviour is present in their sample. Therefore, as far as we are aware, 
these nine galaxies are the only cases where a direct comparison be¬ 
tween the spatial distribution of the blue GCs and the X-ray emis¬ 
sion has been made and published. Excluding NGC 720, all of them 
are massive and giant ellipticals in groups and clusters. 

[Forte et al] d2007Lr2009l . l201^ suggested that the red GCs are 
associated with the bulge stellar component in early type galaxies, 
while the blue GCs would be linked to the halo component. In this 
context, it is natural to think that blue GCs and hot gas can share 
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Figure 18. Azimuthal distribution for the entire sample of GC candidates 
(upper panel), and red, green and blue subpopulations (bottom panel), 
within a ring around the galaxy as described in the text. In order to avoid 
overlapping, the curves corresponding to green and red candidates were 
vertically shifted adding 18 and 35 to the counts. 

the same gravitational potential in equilibrium l lForbes et al.l2012[) . 
However, it is important to notice that similar projected 2-d distri¬ 
butions do not necessarily imply similar 3-d spatial distribution. 

Thus, we decided to make a comparison between the slope 
of the blue GCs density profile and that of the emitting hot gas. 
Fortunately, there is a complete analisys of the X-ray emission 
around NGC 6861 and its neighbour NGC 6868 in the Telescopium 
gal axy group presen ted bv iMachacek et al] ( l2010l) . For both galax¬ 
ies [Machacekrt)^ quantify the X-ray emission through double 13- 
model profiles. The inner component could be associated with the 
galaxy emission and the outer one with a subgroup hot halo. In this 
context [Machacek et al.l suggested that NGC 6868 and NGC 6861 
may each be the dominant galaxy in a merging galaxy subgroup. 
The /3-model is commonly used to re produce the observed X -ray 
surface brightness profiles in galaxies Jo’ Sullivan et al.ir2003l) . and 
is defined as: 

S{r) = S-oll -f (7) 

The relevant parameters in the model are the core radius (rA and 
the slope profile (/3a;), this latter defines the beha viour of the den- 
sity profile for r>rc. In the case of NGC 6861. [Machacek et all 
found the following /3-model parameters: rc=0.5(11.9) kpc and 
/3=0.61(0.38) for the inner (outer) /3-model components, respec¬ 
tively. 

In order to quantitatively compare the slopes of blue GCs and 
X-ray profiles, we have counted blue cluster candidates in concen¬ 
tric circular annuli. Figure [T^ shows counts as small filled circles 
and the corresponding Poissonian error bars. As a red solid line 
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Figure 19. Globular cluster density profile for the blue subpopulation (filled 
black circles) and surface brightness profile in X-ray (red line) ofNGC 6861 
(Jones et al. in prep). The profiles were shifted on the vertical axis for a bet¬ 
ter comparison. Blue dashed line shows the single /J-model fitted to the GCs 
density. The agreement between the blue GCs and X-ray profiles outside 10 
kpc is remarkable. 

we plotted the double /1-model from Jones et al. (in preparation), 
shifted vertically to have a good match. As it can be seen, the X-ray 
profile shows an excellent agreement with the blue globulars den¬ 
sity for r>50 arcsec (> 10 kpc). Subsequently, we have fitted the 
density profile considering points in the range 10.6 <r< 67.1 kpc 
by employing a single /3-model. The obtained val ue of the slope is 
P ~ 0.38 — 0.44, which is very similar to that of iMachacek et al.] 
However, due to the restricted radial range included in the fit, the 
Tc parameter results not well constrained. In order to obtain a more 
reliable P value, we have fitted agai n the data, but k eeping fixed 
the Tc parameter at the value given bv iMachacek et al] We have ob¬ 
tained P — 0.42 ± 0.01. We plotted this fit as a dashed blue line 
in Figure [To] Again, the agreement is excellent showing that outer 
blue GCs present a similar projected distribution to that of the sub¬ 
group X-ray component. 

3.8 Luminosity Function 

The adopted distance of 28.1 Mpc for NGC6861 in this work, 
translates in a distance modulus of (m — M) = 32.24 mag. There- 
fore, taking an approximated tur n over of Mv =— 1.5 mag. (e.g. 
lHarrisl2()0lJ) . and equation 2 from iFaifer et al.l ( 1201 ih . we estimated 
a value for the peak of the GC luminosity function (GCLF) in the 
g' band at -^25.1 mag. As we showed in Subsection 12.41 the com¬ 
pleteness of our photometry is greater than 80 per cent for Qq <25.5 
mag. This indicates that we will be able to obtain a good represen¬ 
tation of the integrated GCLF for NGC 6861. 

In order to do that, we counted candidates in bins of 0.25 mag. 
This raw luminosity distribution is shown in Figure [20] as a dotted 
line histogram. Subsequently we have corrected the counts for com¬ 
pleteness effects and background contamination and plotted it as 
solid histograms in Figure]^ This figure shows that a peak around 
Qq ~ 25 mag is clearly detected. 

Altho ugh different option s exist in the literature to model the 
GCLF (see [Jordan et al.ir2007l) . we have adopted here a Gaussian 



Figure 20. The dotted histogram shows counts of globular cluster candi¬ 
dates considering bins of 0.25 mag. After applying completeness and back¬ 
ground con'ections, we obtain the solid line histogram. The adopted nornial- 
ized background is shown as a dashed histogram. The solid line represents 
the Gaussian fit to the corrected globular cluster luminosity function. 


form. The main reason behind this election is that ther e is a historic 
use of Gaussian represent ation as distance indicator d Jacoby et al.l 
ll992l:IVillegas et al.l2010l) . Therefore, we fit a Gaussian function to 
our histogram, and determined the position of the turn-over magni¬ 
tude (TO), and the dispersion of the GCLF. In the fit we included 
all the bins with magnitude g'<26.2 mag. The values obtained for 
both parameters are listed in Table]^ and Figure[^shows the ob¬ 
tained fit. The value of a'g is in very good agreement with the rela¬ 
tion between the intrinsic Gaussian dispersion of the g-band GCLF, 
and the absolute blue m agnitude of the pare nt galaxy, Og-Ms, pre¬ 
sented in the figure 9 of ilordan et all | |2007|) . 

Our next step was to obtain the GCLF for the different sub¬ 
populations detected in NGC6861. Figure!^ shows the blue and 
red corrected histograms and the fitted Gaussian models. As we 
did not detect any difference between the fit for the green and red 
candidates, we considered both subsamples together in the fit. The 
values of the obtained parameters are listed in Table]^ 

It is interesting to notice that we have obtained a significant 
difference between the turn over magnitudes of blue and red clus¬ 
ters. We can see from Table[^that TOred — TObiue ~ 0.48 mag. 
It is immediate to think in the existen ce of metallicity effects be¬ 
hind this result. [Ashman et al.l (1 1 9951) have investigated the shift 
in the peak of the GCLF as a function of metallicity for the B, 
V, R, I and J filters. They based their results on simulations that 
included the cluster masses drawn from a parent population sim¬ 
ilar to that of the MW, and a mean metallicity ranging from -0.2 
dex to -1.6 dex. According to their table 3, and assuming mean 
\Fe/H] of -1.6 and -0.6 dex for the blue metal poor and red metal 
rich GCs, respectively, we can expect a turn over shift of around 
0.3 mag in the V band. Red GCs in NGC 6861 are probably more 
metal rich than those of the MW. Therefore, the turn over shift in 
this SO galaxy could be larger than in our galaxy. Unfortunately, 
[Ashman et al.l (1 1995!) did not include the Sloan filters in their cal¬ 
culations. For that reason we transformed the t urn over magnitude s 
for blue and red GCs using equation 2 from iFaifer et^ ( 1201 ll) . 
and adopting as mean colours 0.78 and 1.13 mag for blue and red 
candidates, respectively. This translates in a 0.43 mag shift for the 
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a large galactocentric distance. Besides, that value is in agreement 
with seve ral wide field studie s of different G C systems of massive 
galaxies teassino et alj|200^ : lHarriil2009bh . In the inner region 
(<30 arcsec), we assume that the density of GCs is constant. This 
gives us 160 more GCs. 

These 1994 GCs, according to the TO value obtained for the 
whole sample of candidates, correspond to 67 per cent of GCs 
brighter than go=25.5 mag, leaving 33 per cent fainter. Consider¬ 
ing a ±20 per cent variation in the extension of the surface density 
profile due to the uncertainties in their parameters, by integrating 
the LF profile over all magnitudes, we estimate a total number of 
3000 ± 300 GCs. 

Using the apparent visual magnitude U^= 10.97 (NED; an 
identical value was published bv lHo et^l201lh and the distance 
module (m-M)—32.24 for NGC6861, we obtain as its absolute 
magnitude Mv=— 21.21. From this value and from the total num- 
be r of GCs, we can calcula t e the specific frequency {Sn), defined 
bv iHarris & van den Berghl ( ll98lh (Equation [|]l. 


Figure 21. Conected luminosity functions for blue and red subpopulations. 
Dotted lines and vertical short lines represent the Gaussian fit and TO mag¬ 
nitude, respectively. 


Table 8. Parameters obtained by fitting Gaussian functions to the g'-band 
GCLF. Turn-over magnitude (TO) and dispersion (a) for the whole sample, 
and blue and red subpopulations, are shown. 


Population 

TO 

<7 

All 

25.00±0.05 

1.08±0.06 

Blue 

24.70±0.06 

1.00±0.06 

Red 

25.18±0.08 

1.07±0.08 


turn over magnitude in the V band. If we assume that blue GCs in 
N GC6861 are sim ilar to those of the MW, according to figure 7 
of lReiku^ ( |2012|) . the mean metallicity of red GC in this galaxy 
should be ~ -0.25 dex. 

This result shows that it may not be a good idea to measure dis¬ 
tances using the GCLF obtained in blue bands for the whole sam¬ 
ple. It is possible to minimi ze the metallicity ef fect by observing in 
the near IR regime (like g. lVillegas et aljEoiOh . However, another 
option is to use the turn over magn i tude o f the blue GCs as a dis¬ 
tance indicator. IPi Criscienzo et al] (l2006h have shown that adopt¬ 
ing a common EMC fiducial distance scale of (m — M) = 18.50 
mag, the TO magnitudes of metal-po or GCLE in the MW , M31 
and the more distant galaxy sample of lUarsen et al.l (1200 ih are in 
perfect mut ual agreement with Mv=— 7.66 mag. 

Taking lPi Criscienzo et al.l ( l2006h turn over magnitude for the 
blue GCs, the correspond ing mean (p'-i') colours from Table |5] 
and using equation 2 from iFaifer et al. (1201 ll) . we got Mg/=—7.39 
mag. With this value and the apparent magnitude for the turn over 
of the blue GCLF from Table[^ we got a distance modulus of (m — 
M) = 32.1 mag. T his value is in very good agreement with that of 
iTonrv et^ ( l20()lh . (m — M) = 32.24 ± 0.36. mag. 


3.9 Total number of GCs and Specific Frequency 

In Section 13.51 we determined the density profiles for the differ¬ 
ent subpopulations and for the whole system. By extrapolating the 
density profile from ~30 to 600 arcsec, we got 1834 GC. This 
last value is equivalent to ~100 kpc at the adopted distance of 
NGC6861. We have taken this external cut because the analysis 
presented in that subsection showed that there are GCs even at such 


Sn = (8) 

Therefore, for NGC6861 we obtain a surprisingly high value 
of Sn = 9.2 ± 2 . 2 . On the other hand, if we adopt the distance 
modulus obtained from the TO magnitude for blue GCs in previ¬ 
ous section, we get Mv=—21.12 mag, and a Sn = 10.6 ± 2 . 1 . 
These are extremely high values of Sn, simi lar to that of gian t 
cluster dominant early type galaxies like M87. [Harris et"^ j2013h 
presented a compilation of global properties of GC systems. We 
can see that there are only two SO galaxies with D < 100 Mpc 
and Mv<—20 mag which present Sn>1 in that sample. They 
are NGC 1400 and IC3651. However, checking the s ource of that 
values , we can see that in the case of NGC 1400, IPerrett et al.l 
(Il997h found Nt = 922 ± 280 GCs, and a Sn = 5.2 ± 2 . 
On the other hand, iMarln-Franch & Anaricid (l2002h list a value 
of Sn = 5.7 ± 1.8 for IC3651. However, this relies on an indi¬ 
rect estimation of Mv for that galaxy as it does not appear in the 
RC3 catalog, and on an indirect estimation of Nt using the Surface 
Brightness Fluctuation method. In this context, NGC 6861 seems 
to be a more definite case of a high Sn lenticular galaxy. 

High Sn values are usually found in low and high mass early- 
type galaxies. This is evident in the “U” shape plot of S'jv vs My 
like that presented in Fig. 10 in [Harris et alj ( l2013l) . As Sn mea¬ 
sures the formation efficiency of GCs relative to field stars, high 
Sn galaxies can be understood as “field-stars deficient” or “cluster 
rich” systems. However, several studies pr ovided growing evid ence 
in favour of the first situation. Specifically. [McLaughlin! ( 119991) has 
shown that the enhanced Sn value for the most massive galaxies in 
his sample can be accounted for if the mass in GCs is normalized 
to the total baryonic mass of the host (i.e. X-ray emitting gas and 
the stellar component), e = Macs /Mb ~ 0.0026. In this context, 
according its Sn value, NGC 6861 could have a massive baryonic 
halo associated with the hot gas emitting in X-ray. 

Recently, a more fundamental ratio was defined: the ratio of 
the mass in GCs to the total galaxy halo mass, 77 = Mgcs/Mh. 
This ratio seems to be essentially constant over the wh ole mass 
range of galaxies ( [Georgiev et al.l20l^ :[Harris et al.l2 01^. U sing a 
comprehensive database of 307 galaxieZ iHudson et al.l ( [201^ con¬ 
firmed this suggestion through a new calibration of galaxy halo 
masses entirely based on weak lensing. They obtained a value of 
<ri>=3.9 ± 0.9x 10“®. Following Harris et al. (2013) we esti¬ 
mated the total mass enclosed in the GCS of NGC 6861 (adopt¬ 
ing a {M/L)v=2), Mgcs=9 ± 1x10® M©, and therefore we 
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used the < 77 > value to estimate the total mass of the halo of 
NGC6861. We found Mh=23 ± 0.5x10“Mq. This value indi¬ 
cates again that according its GC system, we can expect a massive 
halo in NGC6861. 


4 CONCLUSIONS 

Using deep images taken with GEMINI/GMOS, we obtained multi¬ 
colour photometry of three fields around NGC6861. This galaxy 
presents unique features for an SO type, summarized as follows: 

• NGC6861 presents signs of past interaction or merger, like 
the non-concentric isophotes shown by the galaxy. In addition, we 
have found several different structures that become visible when 
the ratio between the original image and the smooth ellipse model 
is analyzed. 

• The galaxy has a large GC system which extends up to 100 
kpc. By analyzing the colour-colour and colour-magnitude dia¬ 
grams, we detected at least two subpopulations of GC, blue and red, 
with colour peaks at {g'-i')o ~ 0.79 and {g'-i')o ~ 1.07 mag, re¬ 
spectively. By examining the colour histograms and the CMD, we 
inferred the possible existence of a third subpopulation with inter¬ 
mediate colours ({g' — i')o 0.95). We have tentatively separated 
the cluster candidates in blue, green and red subpopulations. How¬ 
ever, spectroscopic and/or IR data are needed to perform an age- 
metallicity study, and confirm or rule out the existence of this third 
subpopulation. 

• The presence of a blue-tilt is detected in the blue subpopula¬ 
tion, probably being the first case in this galaxy type. 

• Regarding the spatial distribution of the blue and green sub¬ 
populations, their density profiles show similar slopes. 

• Outside 10 kpc, the blue subpopulation shows a remarkable 
similarity between the projected density profile and the X-ray sur¬ 
face brightness profile of the galaxy. This result would indicate that 
NGC 6861 co uld be the dominant ga laxy in a galaxy sub-group, as 
mentioned bv iMachacek et al] 1 I 2 OIOI) . 

• The azimuthal distribution of the red and green GC candidates 
exhibit strong signs of elongation, with an ellipticity slightly larger 
than that of the galaxy. 

• By comparing the galaxy light with the red subpopulation, we 
can see that in the inner region they both share similar features such 
as PA, colour and slope of the density profile. 

• Red GC candidates with colours {g'-i')o>l-l6 show an 
asymmetric spatial distribution around the galaxy. Similar features 
were reported in other galaxies, suggesting accretion and/or merger 
of a lower mass neighbour. 

• We obtained the global GCLF of NGC 6861, as well as for 
the blue and red subpopulations separately. A significant difference 
between the TO values of both subpopulations in the g' band, was 
obtained {TOred — TOuue ~ 0.48 mag). The existence of metal- 
licity effects might be behind this result. 

• Using the blue GC candidates, we estimated a new distance 
modulus ( m — M = 32.1 m ag), which shows a good agreement 
with that of iTonrv et alj ( 1200 ih . 

• We have estimated a total population of 3000 ± 300 GCs for 
NGC6861. 

• From the distance given by iTonrv et^ l l200lh . we obtained 
a high specific frequency of Sn ~ 9.2 ± 2 .2, and from our own 
distance estimation, a higher value of Sn = 10.6 ± 2 .1. As far as 
we know, this is the highest Sn measured for a lenticular galaxy 
up to date. This high Sn value could be interpreted as a sign of the 
presence of a massive halo in NGC 6861. 
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